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ABSTRACT 

We present a time-dependent radiative model for the atmosphere of extrasolar planets that takes 
into account the eccentricity of their orbit. In addition to the modulation of stellar irradiation by the 
varying planet-star distance, the pseudo- synchronous rotation of the planets may play a significant 
role. We include both of these time-dependent effects when modeling the planetary thermal structure. 
We investigate the thermal structure, and spectral characteristics for time-dependent stellar heating 
for two highly eccentric planets. Finally, we discuss observational aspects for those planets suitable 
for Spitzer measurements, and investigate the role of the rotation rate. 

Subject headings: planets and satellite: individual (HD 80606b, HD 17156b) — radiative transfer 



1. INTRODUCTION 

Observations of close-in extrasolar planets using 
theSpitzer Space Telescope have permitted the in- 
ference of their atmospheric pro perties via detection 
of their emergent infrared light (iDeming et al.l 120071: 
iKnutson et~alTl2007l : iGrillmair et al.l 120071 : iSwain et all 
120081 ). Most investigations of this type are not sensitive 
to the time variation of the planet's properties. However, 
a real planet will sometimes produce a dramatically time- 
variable infrared signal in respons e to several variables, 
as wa s recently demonstrated by lLangton fc Laughlin I 
(|2008l ). For planets on eccentric orbits, the observed 
signal will vary due to the star-planet distance, as well 
as rapid changes in our viewing angle to the planet. 
The latter is sensitive to the pseudo-synchronous rota- 
tion rate near periastron for planets on eccentric or- 
bits. Also, there is the finite radiative response time, 
wherein the planet's atmosphere does not heat or cool 
instantaneously in response to variable stellar irradia- 
tion. Finally, a vigorous hydrodynamics is expected to 
re-distribute the heating by stellar irradiation. 

A full treatment of time-dependent radiative transfer 
coupled to nonlinear hydrodynamics is beyond the cur- 
rent state of the art. In this paper, we concentrate on 
the effect of time-dependent radiative transfer in isola- 
tion from the re-distribution of heat by dynamics. There 
are two motivations for investigating time-dependent ra- 
diative effects in relative isolation, coupled only to the 
changing rotational aspect geometry. First, we want to 
know the time scales for the exoplanet atmosphere to 
respond as a function of orbit eccentricity and depth 
in the atmosphere. Second, a comparison of our cal- 
culations to Spitzer time series data for eccentric plan- 
ets might in principle inform us of the degree to which 
the observations are affected by dynamics, i.e., can sim- 
ple time-dependent radiative transfer and changing as- 
pect geometry account for most of the observations, or 
do hydrodynamic effects become dominant on a short 
time-scale? Planets on eccentric orbits receive a highly 
varying flux from their stars. Even in the case of one of 
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the mildest ecce ntric planets WASP- 10 b, with an eccen- 
tricity of 0.057 ([Christian et al.l [200 9) this represents a 
variation ~ 25% in the received flux between periastron 
and apastron. Our model takes into account the time- 
variable stellar irradiance impinging on the atmosphere 
of these planets. 

In Section [21 we present our time-dependent model cal- 
culation and assumptions. In Section [3j we will point 
out some aspects of the eccentric exoplanet s. We will 
focus our study on the two transiting planets that have 
the largest eccentricities: HD 80606b (e = 0.932) and 
HD 17156b (e = 0.67). A summary of the results and a 
conclusion will be presented in Section [U 

2. THE ATMOSPHERIC MODEL 

2.1. Radiative transfer 

The energy equation relating the evolution of the tem- 
perature profile in an atmosphere under hydrostatic equi- 
librium is: 

^- = h(p)-c(p), (1) 

where h(p) is the heating rate and c{p) the cooling rate. 
Assuming that the energy flux is purely radiative, we 
have h(p) = -f^^f and c(p) = -f^^ff • We divided 
the net flux into the thermal flux emitted by the atmo- 
sphere Fir (upward — downward) and the net stellar 
flux F* (downward — upward), so that F = Fir — F*. 
The numerical method used to solve the energy equa- 
tion as well as th e sourc es of opacities is discussed in 
detail in llro et al.l ([2005). Once the radiative solution 
is reached, the (temperature-dependent) adiabatic lapse 
rate is calculated for each layer and it is imposed to the 
final solution where it would be exceeded. 

2.2. Input parameters 

Table [1] summarizes the input parameters pertaining 
to each planet, star and orbit among the eccentric tran- 
siting planets. An exception would be the metallicity, 
and we adopted a solar metallicity for all the models in 
this study. In this table, we included HD 209458b, even 
though its eccentricity is consistent with a circular orbit. 
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TABLE 1 

Planetary and stellar input parameters. 





Orbit 


Planet 


Star 


Refs. 


System 


e 


a 


^min 


^max 


P 




Rpl. 


Mass 


-Pspin 


Type 




m 








[AU] 


[AU] 


[AU] 


[d.] 


[°] 


[Rj] 


[Mj] 


[d.] 




[R©] 






HD 80606 


0.93 


0.432 a 


0.029 


0.835 


111.45 


300 


i.i a 


4.18 a 


1.72 


G5 


1.0 


0.43 


(1) 


HD 17156 


0.67 


0.159 


0.052 


0.266 


21.22 


121 


1.23 


3.09 


3.80 


GOV 


1.47 


0.24 


(2;3) 


f HAT-P-2 
\ HD147506 


0.52 


0.068 


0.033 


0.103 


5.63 


190 


0.95 


8.62 


1.87 


F8 


1.42 


0.11 


(4) 


XO-3 


0.260 


0.046 


0.034 


0.058 


3.19 


346 


1.22 


11.79 


2.25 


F5V 


1.38 


-0.18 


(5) 


HAT-P-11 


0.198 


0.053 


0.043 


0.063 


4.89 


355 


0.42 


0.081 


3.94 


K4 


0.75 


0.31 


(6) 


GJ 436 


0.15 


0.085 


0.024 


0.033 


2.64 


351 


0.39 


0.070 


2.32 


M2.5 


0.47 


-0.32 


(7-10) 


WASP-14 


0.095 


0.037 


0.033 


0.041 


2.24 


255 


1.26 


7.725 


2.12 


F5 


1.297 


0.0 b 


(11) 


HAT-P-1 


0.067 c 


0.055 


0.052 


0.059 


3.09 


d 


1.225 


0.524 


3.01 


GOV 


1.135 


0.13 


(12) 


WASP-10 


0.057 


0.037 


0.035 


0.039 


3.52 


157 


1.29 


3.06 


3.45 


K5 


0.784 


0.03 b 


(13) 


WASP-6 


0.054 


0.042 


0.040 


0.044 


3.36 


99 


1.22 


0.50 


3.30 


GOV 


0.87 


-0.2 


(14) 


WASP-12 


0.049 


0.023 


0.022 


0.024 


1.09 


286 


1.79 


1.41 


1.07 


F 


1.57 


0.3 b 


(15) 


HD 209458 


0.014 e 


0.045 


0.044 


0.046 


3.52 


80 e 


1.31 


0.69 


3.52 


GOV 


1.46 


0.01 


(16) 



R eferences. — m f Langton fe Laughlin I (I2008T) : (^IGillon et al.l (I2008T) ; f ^lFischer et al.l (I2007T) f4) lLoeillet et al.) (I2008T) 
; f 5) I Winn eTall (l2QQ8f) ; (6) Bakos et al. (2009) (7) De ming et al.l (120071) ; ( 8) IDemorv et al.l (I2007T) ; (9)lBean et al. (2006) ; 
(10) [M aness et al. (2007) ( ll)lJoshi et al. (2009) ; (12) Johnson et al. (20081) : fl3) IChristian et al.T(l2009T ) : fl4) IGillon et al.l 
(2009) ; (15) Hebb et al. (2009) ; (16) Laughlin et al. (200j). 

Note. — The metallicity is shown only for informational purposes since we used a solar metallicity for our models. The 
minimal and maximal distances (resp. d m in and d max ) as well as the rotation period given by Eq. \S\ (P sp in) are calculated as 
a function of the parameters found in the references. 

a Moutou et al. (2009) analyzed the transit of HD 80606 and found slightly different values with little consequences on our 
study. 

b In fact, it is the global metallicity [M/H]. 

c This value is an upper limit but will serve as guideline in this study. 

inn et al.l ( 2008) m easurements wer e also consistent with a zero eccentricity. 
e The measurement by Laughlin et a l. (2005) are also consistent with e — 0. In that case, the argument of the periastron 
would be irrelevant. 



2.3. Time- dependent calculations 

Since we use a time-marching algorithm to solve for 
Equation [TJ we are able to include the variation of the 
incoming insolation. This variation is due to the non 
constant distance of the planet from its star during the 
eccentric orbit, as well as the rotation of the planet due 
to the pseudo-synchronous sate. This incoming flux can 
be written as a function of time as follows: 



Figure [T] shows the orbit for each eccentric transiting 
planets with the parameters taken from Table [TJ 

Another effect that must be included is that an eccen- 
tric planet is expected to be in pseudo-synchronous rota- 
tion due to the strong tidal forces app l ied to the planet 
during the passage to periastron. iHutl (|1981[ ) calculated 
the remnant spin period as follows: 



Fi nc (t,\) 



(I)) ™vw F ^ 



(1 + 3e 2 
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JL P 6 
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where F| nc is the wavelength-dependent incident stellar 
flux at the top of the atmospheric model, R* is the stel- 
lar radius, d(t) is the time-dependent star-to-planet dis- 
tance, is the planetary longitude of the point of the 
atmosphere that is substellar at t = and F®(\) is the 
flux emerging from the star. At each time step, we calcu- 
late the planet-star distance (d(t)) as well as the rotation 
(#(£)), and the subsequent incident stellar flux. 

A current limitation of our method is that the heating 
rates and cooling rates as well as the abundance profiles 
of the atmospheric constituents are calculated only once 
for the initial conditions. We calculate them at the mean 
distance (semimajor axis) under the condition where the 
incoming stellar flux is averaged over the planet. How- 
ever, the influence of the initial thermal profile chosen as 
well as the heating and cooling rates on the final results 
is only minor, as we will discuss in Section HJ 

3. STUDY OF THE PLANETS 
3.1. Eccentric orbit 



The numerical value of the pseudo-synchronous rotation 
period of the eccentric transiting planets is listed in Ta- 
ble m 

3.2. Individual planets 
3.2.1. HD 17156b 

Discovered by Bar bier i et al.l (j2QQ7h . this long period 
(21.2 days) planet has the second largest eccentricity 
(0.67) among the transiting planets. 

Previous studies — llrwin et "all (j2QQ8h performed a photo- 
metric analysis of the s ystem which gave a pproximately 
the same parameters as lGillon et al.l ()2008l ) summarized 
in Table [T] (except for a slightly larger planetary radius) 
and found no transit timing variation due to the presence 
of additional planets in the system. 

They also presented a hydrodynamic model for this 
planet that predicted theoretical light curves for Spitzer 
observations calculated by using the climate model of 
lLangton fc Laughlinl (|2QQ7D applied to the HD 17156b 
parameters. They predicted the planet to star ratio 
F p /F* in the 8 fim band to vary from ~ 1.7 x 10 -4 to 
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Fig. 1. — Representation of the eccentric transiting planets known so far (circa February 2009) with parameters taken from Tabled The 
center circle symbolizing the star is to scale with respect to the orbit. The arrows indicate the direction to the Earth. The dashed circles 
show the semimajor axis distance (orbit without eccentricity). Note the different scales for each planetary system. 



~5 x 10 -4 in the ~ 30 hr following periastron. The au- 
thors used the iHutl ([1981) formula to obtain P sp in ^ 3.8 
days with HD 17156b parameters. 

Thermal structure — In Figure [2j we can show the tem- 
perature of one latitude during an orbit for several pres- 
sure levels; note the ~ 5.6 rotations during one orbit. 
The low pressure levels are very sensitive to the heating 
variations (the temperature variation of an atmospheric 
parcel is more than 500 K at 1 mbar). Deeper than 1 bar, 
the temperature var iations ar e smal ler than 50 K. As pre- 
viously discussed in llro et all (j2QQ5h . for pressure greater 
than 1 kbar, the thermal structure is mainly driven by 
the internal heat flux originating from the planet interior, 
constrained by planetary evolution models. Moreover, 
due to the longer radiative timescale at larger pressure, 
there is a delay between the heating and the temperature 
increases. 

Figure [3] shows all the thermal structures of one lati- 
tude during an orbit. We can see again that deeper than 
5 bar, the atmosphere is only weakly dependent on the 
incoming heating variation. On the other hand, above 
the 100 mbar level, some of the profiles exhibit a tem- 
perature inversion. This inversion is generated by the 
strong heating by the star when the planet is closer — 
near periastron. It should be noted that no change of 
chemistry is taken into account here so the change of ir- 
radiation is sufficient for causing this inversion. The in- 
version reaches an amplitude of ~ 300 K. This inversion 



has been previously inferred for highly irr a diated planets 
([Burrows et al]l20Q6tlFbrtney et al.ll2QQ8h . iFortnev et all 
(|2008[ ) distinguished two types of planets by the degree 
of incoming radiation causing or not this inversion (pM 
class planets experiencing a very high stellar radiation 
exhibit an inversion whereas cooler pL class planets do 
not). As we show here, the highly eccentric planets are 
then likely to transfer from one class to the other during 
the course of their orbit. 

Figure 2] represents an equatorial cut of the atmosphere 
showing the thermal structure as a function of pressure 
for several phases during an orbit. The rotation leads 
to a moving of the hottest zone with time. As a conse- 
quence of the finite radiative timescale, the hottest zone 
is shifted with respect to the maximum of insolation (tak- 
ing place at the substellar point). Moreover, the eccentric 
orbit generates a global decrease of temperature from the 
first panel to the fourth as well as a uniformization of the 
temperature with a lower day to night contrast. During 
the periastron encounter, the high incoming stellar flux 
leads to a thermal inversion (as shown in panels 1 and 5) 
as noted above. Because of the finite radiative timescale, 
the inversion location is not centered at the substellar 
point but is delayed in the direction of the rotation of 
about 10° (non-constant due to the speed difference be- 
tween the rotation and the heat propagation). 

Flux ratio — Figure [5] shows the planet to star flux ra- 
tio as a function of wavelength. For the coldest profiles, 
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Fig. 2. — Temperature of the planet atmospheres, for a single longitude at selected pressure levels (10 3 mbar, 1 mbar and 1 bar) as a 
function of time. The heating is a factor that is applied to the incident flux with respect to an averaged incoming flux model (F^ nc /F°), as 
defined by Equation [2] We let the calculation make several orbital revolution in order to reach a periodic solution (typically 5 for HD 171 56b 
and 16 for HD80606b). Note the different scales in temperature between the two planets. 




400 600 S00 1000 T200 1400 200 400 600 800 
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Fig. 3. — Temperature structures of the atmospheres for several orbital position (lines) and for each time steps (shaded area). Note the 
different scales. 



L = 0.00088360038 t = 0.25123704 i = 0.50159048 




Fig. 4. — Thermal structure of the planet HD17156b atmosphere for selected times during the orbit (t/P rb = 0; 0.25; 0.5; 0.75 and 1). 
The axis represent the pressure in logarithmic scale. The circles from center out indicate the 10 3 , 1 and 10 -3 bar levels. The substellar 
point direction is indicated by the black line, fixed at a single location to facilitate comparison between frames. The white line indicates 
the direction of the Earth. The time is shown as fractions of one orbit. The area contained in the rectangle is the parcel of atmosphere 
whose temperature evolution is plotted in Figure [2] (in pseudo-synchronous rotation). 
An animation of this figure is available in the online journal. 
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the features — in absorption — are mainly caused by H2O 
and CO. As previous ly discussed in the models including 
a thermal inversion (iHubenv et al.l 12003b iBurrows et al.l 
120061 : iFortnev et al.ll2008[ ). the features are seen in emis- 
sion near periastron due to the inversion of temperature. 

Predictions for Spitzer observations — Croll et al. have 
observed the system during the Cycle 5 of Spitzer with 
IRAC for 30 hr (pid: 50747). The data have not been 
published yet. 

Figure El represents the predicted light curve in the 
ffpitzer 8 /im ban d. It can be compared to Figure 4 of 
llrwin et all (|2008[ ). Our results have a smaller amplitude 
of variation, but we do not account for winds. The flux 
ratio in our model goes from ~ 1.3 x 10 -4 to ~ 4.3 x 10 -4 
during the orbit. 

3.2.2. HD 80606b 

HD 80606b has been discovered bv lNaef et all (|2001f ). 
The high eccentricity of 0.932 leads to a factor ~ 800 of 
insolation between apastron and periastron. 

Previous studies — lLangton fc Laughlin I (|2008f ) per- 
formed a hydrodynamic model of this planet. They pre- 
dicted the planet to star ratio F p /F* in the 8 /im band 
to vary from ~ 4 x 10 -4 to ~ 6.5 x 10 -4 from ~ 20 hr 
preceding periastron unt il ~ 40 hr following p eriastron. 

As recently published, lLaughlin et al . (2009) observed 
the secondary eclipse of this planet with Spitzer in the 
8 /im band. The analysis of the light curve is consis- 
tent with their pre dicted radiativ e time at the 570 mbar 
level of - 4.5hr. iMoutou et all (|2009l ) confirmed that 
the planet transits its star and analyzed a transit event. 
They inferred a planetary radius of R p \ = 0.86 ± 0.1 Rj. 

Thermal structure — The temperature contrasts from day 
side to night side are less pronounced than in the case 
of HD 17156b, due to a globally cooler planet and thus 
a longer radiative timescale as well as a shorter rotation 
period of ~ 1.72 days as we can see in Figure O 

Moreover, as we can see in Figure [3j the planet is glob- 
ally more isothermal below 10 mbar than in the case of 
HD 17156b. The temperature difference from 10 mbar 
to the highest pressures is less than 100 K. On the other 
hand, above 10 mbar, the temperature inversion — when 
present — is very large, reaching a 400 K peak. This is 
due to the fact that although the planet is on average 
cooler than HD 17156b, its periastron distance (~ 0.029 
AU, almost half of HD 17156b's) is so short that during 
the periastron encounter the planet receives a blast of 
radiation from its star. 

Figures [2] and [7] show that the very large planet-star 
distance variation leads to a greater temperature differ- 
ence between periastron and apastron. In Figure El we 
can see that the relatively fast rotation cause a large 
shift between the substellar point and the hottest spot — 
almost 180° in the first panel but highly varying since 
the planet rotates faster than the temperature propa- 
gation and also the radiative timescale depends on the 
atmosphere temperature. 

Our consistent radiative transfer model how- 
ever, gives tempera t ures lower than presented by 
lLangton fc Laughlin I (|2008[ ). We reach temperatures 
up to 900 K, to be compared with the 1200 K obtained 
by these authors. In comparison, our model cannot 



represent latitudinal atmospheric flows but it has a con- 
sistent vertical structure. Moreover, the long radiative 
timescale due to the cold temperatures of the planet 
prevent very large temperature contrasts to occur due 
to radiative effects alone. 

In Figure we can see that the very important temper- 
ature inversion leads to emission features that are very 
large. In contrast, the cold temperatures of the distant 
apastron lead to absorption features that are a lot fainter. 

Predictions for Spitzer observations — In Figure [6j we can 
see that HD80606b's light curve is smooth, due to the fast 
rotation of the planet, but contrary to the hy drodynamic 
simulations of lLangton fc Laughlin I (|2008f ). we do not 
achieve a peak of the flux ratio as strong as 8 x 10 -4 but 
rather 2.1 x 10" 4 . 

This small amplitu de variation does not reproduce 
lLaughlin et al.l (|2009) observed light curve during sec- 
ondary eclipse. In this case, the very short but intense 
heating increase cannot generate the observed planetary 
flux if we account for radiation effects alone and we have 
to take into account hydrodynamic processes. 

4. DISCUSSION 
4.1. Role of the initial inputs to the model 

The main limitation of our model is that we do not vary 
the chemical composition of the atmosphere in lock step 
with changes in the thermal profile. In order to quan- 
tify the error produced by using a constant abundance 
profile of the constituents, we conducted a test for the 
planet HD 17156b assuming the two different extreme 
values for heating and cooling rates, as well as the abun- 
dances, calculated for the hottest and coldest profiles ob- 
tained by the previous study. Although the calculations 
take a different amount of time before reaching a peri- 
odic solution (the longer radiative timescale for a colder 
thermal profile makes the model slower to converge), the 
final periodic solution gives similar results for a complete 
orbit. Figure [8] shows the relative differences in the fi- 
nal thermal structure as a function of time using the two 
extremes initial temperature and heating/cooling rates. 
The largest value is less than 0.5%. 

This is consistent with the estimation of the calculation 
time of CO-CH4 chemical reaction, which is the one de- 
termining the main spectral features of the planet. We 
calculated th i s rea ction time using the prescription of 
iBezard et all (|2002f ). For HD17156b, this time above the 
250 mbar pressure level — where the temperature vari- 
ations begin to be relatively important — is less than 
10 T s, more than the time of the orbit. And the cooler 
planet HD80606b has even longer chemical timescales 
^ 3 x 10 15 sec. These results justify our approach of 
not changing the chemical composition with time. 

In order to confirm this conclusion, we calculated the 
flux viewed from Earth integrated in the 8 /im band for 
the two extremal initial profiles for HD17156b. This re- 
sults in little differences, as shown in Figure [9l 

We have adopted solar composition for all our calcu- 
lations. If we increase the metallicity, we expect the 
general thermal structure to be slightly hotter due to 
the highe r opacity and thus g reater absorption of incom- 
ing flux ( Fortnev et al.l [200 6) . Global spectral features 
have been shown to exhibit only minor dif ferences due to 
metallicity in the 8 /im spectral window ([Burrows et al.l 
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Fig. 5. — Ratio of the planet flux by the star flux as a function of wavelength for every times (shaded area) and at several positions 
(lines). 

A color version of this figure is available in the online journal. 
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Fig. 6. — Evolution of the flux ratio ( P F an ' ) in the Spitzer 8 /im as a function of time. We integrated the flux emited from the hemisphere 
seen from Earth; the viewing geometry from the orbit as well as due to the pseudo-synchronous rotation has been taken into account. Note 
that the scale changes between the two panels. 
A color version of this figure is available in the online journal. 
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Fig. 7. — Same as Figure [4] for the planet HD80606b (with a different temperature scale). 
An animation of this figure is available in the online journal. 
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120061 : iFortnev et all 12010( 1. On the other hand, the ra- 
diati ye time scale is propor tional to the atmospheric den- 
sity (| Showman et all 12008( 1 which increases with metal- 
licity, but this timescale also decreases with temperature 
(higher metallicity atmospheres are hotter). These two 
effects should compensate each other to a large extent. 
Consequently, we expect our qualitative conclusions to 
be insensitive to metallicity variations, to first order. 

4.2. Role of the rotation rate 

It should be noted that the formalism of iHutl (|1981l ) 
contains uncertainties in the calcu lation of the rotation 
period of the planet. For instance, llvanov fc Papaloizoul 
(|2007[ ) using another formalism gave a rotation rate equal 
to ~ 1.5 times the circular orbit angular velocity at peri- 

astron Q, p given by: £} p ~ . / ^Ir^S where M* is the mass 

V min 

of the parent star. With the parameters of HD 17156b, 
this gives a rotation period of P^ n ~ 2.5 days. 

In order to test the effects of the period of rotation, we 
used the nominal model for HD 17156b with twice the 
rotation rate and half the rotation rate given by Equa- 
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Fig. 8. — Relative difference of resulting temperature at each 
level between the two extremal initial profiles. Profiles using the 
initial conditions from the hottest and coldest profiles achieved 
by the nominal calculations for HD17156b have been calculated. 
Here, the relative differences between the two extreme cases for 
each pressure level as a function of time are shown. 
A color version of this figure is available in the online journal. 
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Fig. 9. — Relative difference of resulting flux ratio integrated in 
the 8 jum band. As in Figure [8] this shows that using the coldest 
or the hottest profile from the nominal calculations for HD 171 56b 
as initial starting point has little effect on the final light curve 
obtained. 



2 x Pgpfn ~ 7.6 days) as well as the rotation rate given by 
llvanov fc Papaloizoul ([20071 ) (P spin = P^ in ~ 2.5 days). 

Moreover, a more rapid rotation rate can mimic a 
planet wherein uniform zonal winds occur in addition 
to the rotating frame. 

Figures [10] and [11] show the temperature structure of 
the planet (equatorial cuts) as a function of time for the 
two opposite cases. The planet with a rotation rate 2 
times faster than our nominal profile is more uniform, as 
expected. The first and last panels of Figure [IT] exhibit 
temperature contrasts up to 300 K at 1 mbar and as 
the other panels show there is almost no latitudinal con- 
trast below this pressure level when as the planet to star 
distance increases. On the other hand, with the slow- 
est rotation rate as shown in Figure [lOj the day /night 
contrasts are more pronounced, with a 500 K day/night 
contrast at the 1 mbar level near periastron and still a 
~ 50 K contrast far from periastron. 

Figure [12] shows the Spitzer 8/im light curve for 
the three cases. The differences are large enough 
that observations should be able to discrimina te be- 
tween the rotation rate pre dicted by IHutl (jl98l[ ) or by 
llvanov fc Papaloizoul (|2007[ ). 

5. CONCLUSIONS 

We presented a time-dependent radiative model for the 
atmosphere of the transiting planets. This model takes 
into account the eccentricity of their orbit as well as the 
pseudo-synchronous rotation. 

We applied our model to the planets HD 80606b and 
HD 17156b. The results define the radiative portion of 
the response of these eccentric planets to time- variable 
stellar insolation and different possible values of the their 
pseudo-synchronous rotation rate. 

We showed that the high eccentricity makes the planets 
periodically transfer from a state where a stratospheric 
temperature inversion takes place to a cooler planet with 
monotonic temperature profiles and vice versa. Signifi- 
cant time variation in an eccentric planet's spectral fea- 
tures would be expected based on our calculations, and 
methods to observe such variations phased to the orbit 
would be very valuable. 

Moreover, the maximum of temperature and observ- 
able flux is delayed with respect to the maximum of stel- 
lar heating (periastrion) , due to the finite radiative time 

In some cases, an analysis of the planet during the pe- 
riastron encounter is feasible. It will in principle allow 
us to discriminate which of the competing theory that 
predict the pseudo-synchronous rotation period is cor- 
rect since the temporal variation of the planet flux is 
relatively highly dependent on its rotation period. 

Before reaching a time-dependent radiative model in- 
cluding hydrodynamics, such observations will also help 
determine if the planet structure is mainly driven by ra- 
diation or by atmospheric dynamics. 

In the case of HD80606b, the heating increase takes 
place for a such short time that the radiative model fails 
to explain the observed flux variation without accounting 
for hydrodynamic effects. 
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the NASA Postdoctoral Program at the Goddard Space 
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Fig. 10. — Same as Figure [4] for P sp i n = 2 x P^f* This figure is also available as an mpeg animation in the electronic. 
An animation of this figure is available in the online journal. 
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An animation of this figure is available in the online journal. 
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